** 



On Bifurcation Points of the Stationary 
VI asov- Maxwell System with Bifurcation Direction 

o" 

o: N.A.Sidorov W, A.V.Sinitsyn 

<D ' ( \ 

C/0 ; {*> Department of Mathematics 

Irkutsk State University 
Karl-Marks Str., 1 
(N ; Irkutsk 664003, Russia, 

O ' e-mail: sid@math.isu.runnet.ru 

G\ 
O 
O 
O 

-i— > 

B 

> 

Let us consider a multicomponent plasma containing electrons and positively charged 

ions of various kinds and described by the multiparticle distribution function /, = /j(r, v), 
i = 1, . . . , N. Plasma occupies a domain D G R 3 with smooth boundary. The particles 
interact via a self-consistent force fields; collisions of particles are neglected. 

The plasma is described by the following version of a nonrelativistic stationary Vlasov- 
Maxwell(VM) system [1] 

v ■ 8 r fi + qi/mi(E + -v x B) ■ d v f { = 
c 
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r eD GR% v eR 6 , i = l,...,N, 



'' v eR 3 , 
curlE = 0, 
divB = 0, (7) 



N 

A 

v = p, 

R 3 



divE = 47rJ2<lk / fk(r, v)d 
fc=i JR3 

4vr N r 
curlB = — ^qk I vf k {r, v)d 



v = j. 



Here p(r), j(r) are the densities of the charge and the current, respectively, and E(r), 
B(r) are the electric and the magnetic field. 

We shall search the solution E, B, f of VM system (I) for r E D C R 3 with the 
boundary conditions on the potentials and the densities 

U\dD=u i, (A,d)\ dD =Uo2, (2) 

P \dD= 0, j \ dD = 0, (3) 

where E = —d r U, B = curl A, and U, A are the scalar and the vector potentials. 

The solution E°, B°, for which p° = 0, j° = in a domain D, will be referred to as 
the trivial one. 

Our aim is to obtain the existence theorem of nontrivial solutions for the boundary- 
value problem (I), (2), (3). 

In this paper we restrict ourselves to distributions of the form 

/•(A, r, v) = a{\)U-<X)v 2 + <Pi(X, r),v d t (X) + ^(A, r)) = o(A)/,(A, R, G), (4) 

<Pi : R 3 > R : tpi : R 3 ► R, r E D C R 3 , v E R 3 ; 

\ER, aER + , Oil E R + = [0,oo),; d { ER 3 , i = l,...,N, 

where the functions tpi, ipi generating the corresponding electromagnetic field (E, B) are 
to be determined. Let us note that similar distribution functions were introduced in [2] 
and have been used in [3]. We are interested in the dependence of desired functions tpi, 
ipi generating the nontrivial solutions on the parameter A in distribution (4). The case, 
when the parameters a iy d t are independent of A was considered in [4-6]. Here we consider 
the general case. 

Let us give a one preliminary result on reduction of the VM system (I) with conditions 
(2) to the qusilinear system of elliptic equations for distribution (4) [7]. We assume that 
the following condition is satisfied: 

A. /j(R, G) in (4) are fixed, differentiable functions; ««, di are free parameters, | di |^ 
0; ifi = Cu + h(p(r), tpi = c 2i + kiip(r), cu, c 2 i are constants; the parameters /j, hi are 
related by equations 



the integrals 



converge for all tpi, ipi 



a±qi rrii m\ rrii 



fidv, / fiV dv 

R3 JR3 



Theorem 1. Suppose that the distribution function f\ has the form (4) and condition 
A is satisfied. Let the vector-function (tp, ip) be a solution of system of equations 



^naq 



Ay? = a(A)/i Y^ Ik \ fk(X)dv, fj, 

i. i J R3 



i i J&~ ' rn 



4vrg 
if, v — - 



Aip = a(\)v J2<lk (v, d)f k (\)di 

i.i J R? 



,. , Jr? " mc 2 



J = — y^Qk vfkdv, (3 - const. 



To this solution there correspond potentials 



with conditions (2). 
The proof see in [6] . 
We introduce the notations 



ji = / vfi dv, pi= fidv, i = l,...,N. 
Jr 3 Jrz 

and impose the following condition: 

B. There exist vectors fa G R 3 such that ji = fapi, % — 1, . . . , N. 
Let condition B holds, then system (5) becomes 

N 

Ay? = a(\)aJ2QiA, 

i=l 



(5) 



<P \dD= Uqi, ip \ dD = — u 02 (6) 

m mc 

in the subspace 

(d r <pi,di)=0, (d r ij; i ,d i )=0, i = l,...,N. (7) 

Then the VM system has the solution 

m d /"l mc 

E= 2^q dr ^ B = ^(P + J (dxJ(tr), r )dt)-[dxd r ^ — , (8) 

where 
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u = -■£?>> A=^i^ + A 1 (r) 7 (A u d) = (9) 



(10) 



N 

Aip = a(X)u^2q i (fa,d)A i , 
i=i 

where 

MKh<P,kitl>)= I hdv. (11) 

From now on, for simplicity, we consider the auxiliary vector d directed along the axis 
Z. Then ip = ip(x,y), ip = ip(x,y), x, y G D C R 2 in system (10). Moreover, let iV > 3 
and j 1 7^ const. 



Let Dbea bounded domain in R 2 with boundary dD of class C 2,a , a G (0, 1]. 
We introduce a continuous vector-function of parameter A 

e(A) = {kip , M°, «i, di, . . . , l N tp°, k N ^, a N , d N ) e R 4N (12) 

and a contraction of the set (the value of this vector) induced by the vector (12) and the 
boundary conditions (3) for the local densities of the charge and the current 

EfeLi qkA k (l k (p°, k k i/>°, Oi k , d k ) = 

n =M „ } (13) 

E fe =i Qk(Pk, d)A k (l k (p°, k k ip°, a k , d k ) = C 



with 



(p o = _^g ^o = J_ Mq2? N > 3 (14) 

m mc 



We introduce the condition. 

C. Let e(X) G Q for VA from the some open set of real axis for x G (— r, r). 
Then the system (10) with the boundary conditions 

<P \dD= <P°, 4> \dD= ^° (15) 

has the trivial solution <p = tp°, if) = ip° for VA G (—r,r) by (13), (14). By Theorem 1, the 
boundary- value problem (1), (2), (3) has the trivial solution 

77? 

E° = d r <p° = 0, reDcR 2 , B° = (3d 1 , 

/° = a(A)/,(-a,(A)7; 2 + c u + h<p°(\), (v, d,(\) + c 2l + k^°(X)). 
for all A G (— r, r). Moreover, p° = and j° = in a domain D. 

Definition 1. A point e = e(X ) G fl is called a bifurcation point of solution of the 
VM system (I), (2), (3) with a bifurcation direction e = e(X), where A : (—r,r) — > Q be 
the continuous vector-function of A, if every neighborhood of the vector (e°, E°, B°, f°) 
corresponding to the trivial solution with p° = 0, j° = in a domain D contains a vector 
(e*,E,B,f) with e* = s*(X), X G (-r,r) satisfying system (I) with conditions (2), (3) 
such that 

|| E - E° || + || B - B° || + || / - f° ||> 0. 

Moreover, the densities p and j interior to domain D need not vanish. 

Using the Taylor expansion and singling out the linear terms, we rewrite system (10) 
in operator form as follows: 

[L - o(A)Li(e(A))]u - R(e, u) = 0, (16) 



where e = e(X) be a bifurcation direction of (12). Here 



A 
A 



u=(<p-<p°), V-V' ) 7 



A? 



L i = J2 1* 



s=l 



III ^ 

/■*'« dx 



t^ks-gf 



dA s 



vl s ((3 s ,d)^ vk s ((3 s ,d)^ 



A 



x=l s if°, y=k s tp° 



z/T 3 vT^ 



Operator R(e(X),u) is analytic in neighborhood of the point u — 0: 



R(e,u) = J2J2 e™( u ) b s 



i>l s=l 



where 



A q. 



d 







n ' dx 
are ith-order homogeneous forms in u and 

Qil+12 



Qis{u) = -^-(huiTT + k s u 2T -) l A s (l s ip° ) k s tp° 



dy' 



dx h dy i 



A„(e, x, y) \ x= i sip o } y=ks ^o= for 2 < i x + i 2 < I - 1, s = 1 



(//, z/(/5„rf))/. 



(17) 



,...,JV, 



The existence problems for a bifurcation point, e G fi, of the boundary-value problem 
(10), (15) can be restated as the existence problem for a bifurcation point for the operator 
equation (16). 

Let us introduce the Banach spaces C 2,a {D) and C°' a (D) with the norms || • ||2 jQ , and 
|| • ||o, a , respectively, and let W 2 ' 2 (D) be the ordinary Sobolev L 2 - space in D. 

Let us introduce the Banach space E of vectors u = (■u 1 ,-u 2 ) 7 ! where Ui G L 2 (D); L 2 
is the real Hilbert space with inner product ( , ) and the corresponding norm|| • ||l 2 (d)- 

A ° 

We define the domain D(Lq) as the set of vectors u = (ui,u 2 )'with U{ EW 2 ' 2 (D). Here 

o 

W 2 ' 2 (D) consists of W 2,2 functions with zero trace on dD. 

Then L : D C E — > E is a linear self-adjoint operator. By virtue of the embedding 

W 2 ' 2 {D) C C°' a (D), < a < 1, (18) 

the operator R : W 2,2 C E — > E is analytic in a neighborhood of the origin. By the 
embedding (18), any solution of Eq. (16) in D(L ) is a Holder function. Moreover, 
since the coefficients of system (16) are constant, the vector R(s,u) is analytic, and 
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3D G C 2 ' Q ; it follows from well-known results of regularity theory for weak solutions of 

o 

elliptic equations [8] that the generalized solutions of Eq. (16) in W 2,2 (D) actually belong 
to C 2 ' a (D). 

The operator L\ e L(E — > E) is bounded and linear. Under the above assumptions 
on Lq and Li, all singular points of the operator 

L(e) = L - o(A)L 1 (e(A)) 

are Fredholm points in sense [9] . 
Let us introduce the conditions 

I. 7\ < 0; 

II. T X T A - T 2 T 3 > 0. 

If -rp- \ x=y=0 > 0, then inequality II is satisfied. Let us introduce the matrix 
II Qij \\i,j=i,...,N, where 0^ = qiq^ljh - kjk)^ - (3 h d). If the derivatives ^, ^ are 
positive and equal at the point x = y = 0, 0^- > 0, % ^ j, then conditions I, II are 
satisfied. Evidently, the elements of 0^ are nonnegative, because of identities 



N. 



q% (di, d) d 2 hi 
sign— = signq, = — — , q < 0, g« > 0, % — 2, 

Let us denote by E the matrix generating the operator Li. 

Lemma 1. Let conditions I, II be satisfied. Then the matrix E in (17) has two simple 

eigenvalues 

87TQ; | q | 



X+(e) =//Ti + o(l),// 



/;/ 



X-{ £ ) =V 



T ± T 4 - T 2 T 3 1 



o e 



r/ 



47r | g 



>0 



Ti <r m 

for 4 -»• 0. 

The eigenvectors corresponding to negative eigenvalue X- °f the matrices S and E', re- 
spectively, are 



("2 



2k 



+ 0( 



1 





1 



Ol 



(19) 



Proof see in [6, Lemma 4]. 

Let /xo is a n-multiple eigenvalue of Dirichlet problem 

-Ae = /xe, e \ dD = 0, 



(20) 



Let us introduce the condition: 



D. Suppose that condition C is satisfied and there exists Ao G (—r,r) such that 
a(AoX-(e(Ao)) +^0 = 0; moreover, a(A)x_(e(A)) is the monotone increasing (decreasing) 
function in neighborhood of Ao- 

Rewrite the system (16) in the form 

Bu-B 1 (\)u = R(e(\),u), (21) 

where we introduce the following notations 

B = L + fi , Bi(A)=Li(e(A))+ju . 

The operator B is Fredholm and it has n eigenfunctions of the form 

\ ei, i = l,...,n. 
c 2 J 

Vector (ci, c-i)l = c is defined by (19) up to the constant multiplier and Lic=x_c. There- 
fore, 

BxGi = a(X)L 1 (e(Xo))ce i - a(X )L 1 (e )ce i = (a(A)x_(e(A)) - a(A X-( £ o))ce; 

and by taking into account condition D we obtain the following identities 

< Biei,ej >= (a(A)x-(e(A)) + yu ) | c | 2 %, i,j = l,...,n. (22) 

Theorem 2. Let conditions A, B ; C, D ; I, II, as well one of the following conditions 
be satisfied: 

l)/i is odd-multiple eigenvalue of problem (20); 

2) fi = a(A)/j(— oii(v 2 + (fii + v • di(X) + ipi). Then e = e(\ ) is a bifurcation point of the 
boundary problem (1), (2), (3) with bifurcation direction e = e(X). 

Proof. In order that to prove Theorem 2, we need to show that Ao is a bifurcation 
point of Eq.(21). The corresponding finite-dimensional branching system was given in the 
paper [6, Eq.(5.1)] for the case, when /j, dj are independent of A and o(A) = A G R + . It 
is shown in the same place that Ao is a bifurcation point if and only if Ao is a bifurcation 
point of branching system. 

In general case, when the values «j(A), di = di(X) are dependent of A, we can reduce 
the bifurcating system to the form 

(a(X)x-(e(X))+no)(c 2 1 + cl)Z i + r i (Z 1 ,...,Z n ,X) = 0, i = l,...,n. (23) 

with | r(£, A) |= o(\ £ |), using identities (22) and results from [6]. Here the vector-function 
r(£, A) in equation (23) is potential of £, if the condition 2) of Theorem 2 be satisfied. 

Because the continuous function a(A)x_(£(A)) + /xo is equal to zero at the point A , 
and be monotonic in it neighborhood, then in any neighborhood of point Ao, £o — exists 
a couple A*, ^* with £* ^ satisfying the system (23) [ 6, see the proof of Lemma 8 and 
terminology in 10]. Hence, Ao is a bifurcation point of the system (23). Theorem 2 is 
proved. 
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